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Abstract 


The social cost of carbon dioxide (SC-CO;) measures the monetized value of the damages to 
society caused by an incremental metric tonne of CO» emissions and is a key metric informing 
climate policy. Used by governments and other decision-makers in benefit-cost analysis for over 
a decade, 5С-СО; estimates draw on climate science, economics, demography, and other 
disciplines. However, a 2017 report by the US National Academies of Sciences, Engineering, and 
Medicine* (NASEM) highlighted that current SC-CO2 estimates no longer reflect the latest 
research. The report provided a series of recommendations for improving the scientific basis, 
transparency, and uncertainty characterization of SC-CO2 estimates. Here we show that 
improved probabilistic socioeconomic projections, climate models, damage functions, and 
discounting methods that collectively reflect theoretically consistent valuation of risk, 


substantially increase estimates of the SC-CO2. Our preferred mean SC-CO» estimate is $185 per 
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tonne of СО; ($44-413/t-CO:: 5-95% range, 2020 US dollars) at a near-term risk-free discount rate 
of 2 percent, a value 3.6-times higher than the US government's current value of $51/t-CO2. Our 
estimates incorporate updated scientific understanding throughout all components of SC-CO2 
estimation in the new open-source GIVE model, in a manner fully responsive to the near-term 
NASEM recommendations. Our higher SC-CO; values, compared to estimates currently used in 
policy evaluation, substantially increase the estimated benefits of greenhouse gas mitigation and 


thereby increase the expected net benefits of more stringent climate policies. 


Main 

Policies to mitigate greenhouse gas emissions are often evaluated in terms of their net benefits 
to society. The net benefit of a climate policy is the difference between the economic cost of the 
emission reduction (the mitigation costs), and the value of the damages that are prevented by 
that emission reduction (climate benefits, among others). In regulatory impact analysis the 
climate benefits of СО; emission reductions are typically computed by multiplying the change in 
CO» emissions caused by the policy with an estimate of the SC-CO2. This makes the SC-CO» a 


highly influential metric, informing analysis of a wide range of climate policies worldwide. 


For more than a decade, the US government has used the SC-CO» to measure the benefits of 
reducing carbon dioxide emissions in its required regulatory analysis of more than 60 finalized, 
economically significant regulations, including standards for appliance energy efficiency and 
vehicle and power plant emissions?. In the United States, the SC-CO; has also been used as the 
basis for federal tax credits for carbon capture and storage; proposed federal carbon tax 
legislation; state-level zero emission credit payments for nuclear generators and power sector 
planning; among other applications?. The SC-CO2 also supports decision making by government 
environmental agencies in other countries (e.g., Germany, Canada, and Mexico), and is used in 


standardized corporate environmental and sustainability accounting’. 


The SC-CO; is estimated using integrated assessment models (IAMs) that couple together 
simplified representations of the climate system and global economy to estimate the economic 
effects of an incremental pulse of СО; emissions. These models generally follow a four-step 


process in which (1) projections of population and GDP inform a СО; emissions pathway; (2) the 
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СО; emissions path drives а climate model that projects atmospheric greenhouse gas 
concentrations, temperature changes, and other physical variables such as sea level rise; (3) the 
resulting climate change impacts are monetized and aggregated as economic damages; and (4) 


economic discounting combines all future damages into a single present value. 


In 2017, а NASEM report assessing the SC-CO2 estimation methodology used by the US federal 
government found that the leading IAMs used for estimating the SC-CO2 have not kept pace with 
recent advances in climate, economic, and demographic science. The NASEM report offered 
near-term recommendations for improving each step of the SC-CO» estimation process to 
improve the scientific basis, characterization of uncertainty, and transparency of the SC-CO;. 
Recently, Executive Order 13990 re-established the US Interagency. Working Group on the Social 
Cost of Greenhouse Gases (IWG) to update the federal government's official 5С-СО2 estimates, 
and to consider these NASEM recommendations in the process. Others have also criticized the 
models supporting the past federal SC-CO» estimates for a number of problems including 
damages representations that do not reflect recent science, outdated climate system models, 


and imperfect characterization of the compounding uncertainties affecting SC-CO2 estimates." 


Here, we provide probabilistic SC-CO2 estimates from the Greenhouse Gas Impact Value 
Estimator (GIVE), a newly created integrated assessment model designed for quantifying the 
benefits of emission reductions. The model is built on the Mimi.jl platform, an open-source 
package for constructing modular integrated assessment models?. By using novel components 
for each step of the SC-CO? estimation process, GIVE incorporates recent scientific advances that 
аге unaccounted. for by the previous generation of IAMs used in regulatory analysis. Crucially, 
GIVE quantifies uncertainties in each component and propagates these compounding 
uncertainties through the entire computation, thus allowing for a theoretically consistent 


valuation of the risk associated with a marginal emission of CO2. 


Each individual component in GIVE is based on recent peer-reviewed research on socioeconomic 
projections, climate modelling, climate impact assessments, and economic discounting. We 
implement GIVE with a set of internally consistent, probabilistic projections of population?, per- 
capita economic growth??, and CO», CHa, and N20 emissions? generated using a combination of 
statistical modelling and expert elicitation, collectively referred to as the Resources for the Future 
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Socioeconomic Projections? (RFF-SPs). Many existing |АМ5 use outdated climate models and have 
been shown to produce temperature dynamics inconsistent with more sophisticated Earth 
system models*!!, Further, damage functions supporting previous SC-CO2 estimates аге, to a 
large extent, based on studies from several decades ago!. A vast literature since then has 
expanded and improved our scientific understanding of how changes in climate will likely affect 
human wellbeing". To address these shortcomings, we combine socioeconomic uncertainty with 
probabilistic models for the climate system and damage functions (defined as functions that 
relate changes in climate outcomes such as temperature to economic impacts in dollars). The 
GIVE model employs the FalR v1.6.2 climate тоде! 33, the BRICK sea-level model? 7, and 
updated damage function components representing the latest empirical research for the impacts 


of climate on agriculture!?, mortality’®, energy consumption?9, and sea-level rise?!. 


Recent important contributions to the SC-CO; literature have generated improvements to various 
components used by integrated assessment models??? (see Supplemental Information section 
SI.3 for an overview of this literature). The GIVE model's key contribution to this literature is the 
holistic implementation of recent advances in probabilistic socioeconomics accounting for policy 
uncertainty, fully quantified scientific uncertainty including climate tail risk and sea-level rise, 
addition of non-market sectoral damages (i.e., costs not included in GDP accounting such as 
mortality risk), and economic discounting tied to uncertain economic growth. These advances 
allow for a full valuation of the risk resulting from those compounding uncertainties based on 
improved scientific, economic, and demographic evidence? which have previously been 
unavailable. The GIVE model's implementation of this comprehensive set of scientific 
improvements affirms a key result from recent work on the SC-CO27^?7, namely that improved 
scientific understanding of the components of SC-CO; calculation leads to a higher SC-CO; than 
has been previously used in US policymaking; moreover, our approach demonstrates this using a 
more robust methodology reflecting the current state of the literature. GIVE's inputs and outputs 
are spatially resolved at the level of 184 countries for population, income, and damages (except 
for agriculture damage outputs which are resolved at 16 regions). Climate change has the 
potential to exacerbate existing economic inequities??97?, and our work would allow future 


consideration of this issue through equity weighting??. 
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We calculate the 5С-СО; as the discounted sum of additional damages per incremental tonne of 
CO» produced by an emissions pulse in 2020 along an uncertain emissions trajectory derived via 
formal expert elicitation that reflects continued technology and policy evolution. We use an 
empirically calibrated stochastic discounting framework consistent with the observed behaviour 
of interest rates and economic growth?!. We provide 10,000 SC-CO; values using a Monte Carlo 
approach that samples interrelated socioeconomic, climate, and damage function uncertainties 
(Extended Data Table 2). The GIVE model can also be used to compute the social cost of other 


greenhouse gases (e.g., CHa, N20, HFCs). 


We illustrate the relative importance of our updated model components by comparing them to 
outputs from the well-known DICE model?^. We also assess the sensitivity of our SC-CO> estimates 
to our choice of sectoral, regionally disaggregated damage functions by comparing them to two 
aggregate, global damage functions based on meta-analyses of the broader damages 


literature?2?3, 


Socioeconomic projections of economic growth, population, and greenhouse gas emissions 
represent important sources of uncertainty in the SC-CO». In previous models, this uncertainty 
has been poorly characterized?3*?°, Population and growth scenarios based upon the Shared 
Socioeconomic Pathway (SSP)?9 narratives, which were prominently featured in the 
Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report (AR6)!^, do not 
typically come with associated probabilities, though there have been efforts to assign such 
probabilities a posteriori based on expert surveys??. The small number of SSPs precludes sampling 
the large and continuous space of possibilities that characterizes future socioeconomics and 
emissions. A strength of scenario-based analysis is in the qualitative exploration of uncertainty, 
for example through the use of bounding scenarios, including scenarios accounting for outcomes 
ме outside the range of historical experience that become increasingly possible over very long 
time horizons. Such an approach does not, however, facilitate the quantitative evaluation of 
uncertainty and the calculation of expected values, a common requirement for policy analysis. In 
some cases, a lack of quantification of relative probabilities can lead to disagreements over what 
scenarios constitute a plausible reference саѕе?8 ^9, A holistic, probabilistic approach to 


accounting for these uncertainties was recently іпігоаисеа ^2, Building on this approach, we 
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sample the RFF-SPs, comprising multi-century probabilistic projections of population? and GDP 
per capitat? at the country level as well as a distribution of projections of global CO2, CHa, and 


N20 emissions derived from a combination of statistical and expert-based approaches. 


The RFF-SPs complement the scenario-based approach by providing an alternative approach that 
characterizes the joint uncertainty across annual GDP, population, and greenhouse gas emissions 
for the multi-century timespan required for climate damage estimation. They also leverage 
expert knowledge to account for potential future changes in policy and technology. They project 
that (Fig. 1): median world population peaks at 11 billion around 2130 and subsequently declines 
to 7.3 billion in 2300, (2.8 billion—21 billion: 5-9596 range); median global per capita annualized 
economic growth declines slowly to reach a cumulative time-average rate of 0.8896 between 
2020 and 2300 (0.17%-2.7%: 5-95% range); median net global CO? emissions decline to roughly 
40 percent of today's levels in 2100 (-20%-150%: 5-95% range), with slower declines thereafter 


(see Supplemental Information section SI.1 for more detail on the RFF-SPs). 


Our mean 5С-СО; estimate using the preferred discounting scheme is $185/t-CO2 ($44-413/t- 
CO»: 5-95% range, in 2020 US dollars, as are all dollar results in this study) (Fig. 2). This is 3.6 times 
greater than the US government's current, most commonly cited mean value of $51/t-CO; using 
a 396 constant discount rate?. We report mean 5С-СО; values throughout this paper to align our 
results with the standard expected. net benefit framework that is routinely used for policy 


analysis^ and supported by standard economic theory^?^6, 


SC-CO: estimates are well-known to be highly sensitive to the discount rate?? because the long 
residence time of СО; in the atmosphere means a СО; emissions pulse continues to cause 
damages long after it was emitted. Our preferred discounting scheme uses a 296 near-term risk- 
free discount rate, which reflects the recent literature on real interest rates^" ^?, which have 
declined substantially over recent decades???!, as well as the central tendency from a survey of 
academic economists**. Our discount rate is related to stochastic consumption growth іп a 
Ramsey-like equation, which is the commonly used approach to value marginal impacts amid 
uncertainty in future payoffs and consumption levels?95^, In this way, the parameterization of the 
discount rate captures risk preferences using the risk aversion parameters discussed in the 


Methods. 
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We also assess (Extended Data Fig. 1 and Table 1) the sensitivity of our SC-CO» estimates to 
discounting by also using near-term rates of 396 ($80/t-CO» mean: $12-197/t-CO2 5-95% range), 
to facilitate comparison with the US government's current, most commonly cited $51/t-CO; 
figure, as well as 2.5% ($118/t-CO2 mean: $23-280/t-CO» 5-95% range) and 1.596 ($308/t-CO; 
mean: $94-626/t-CO2 5-95% range). We additionally show (Extended Data Fig. 2) the temporal 
evolution of the discounted marginal damages by year based upon the preferred 296 near-term 


discount rate case. 


Our SC-CO» estimates are based on regionally disaggregated damage functions for four sectors. 
As a sensitivity analysis, we replace the sectoral damage functions in GIVE with two distinct, 
globally aggregated damage functions that are based on meta-analyses of the climate impacts 
literature???3, Under a 2% near-term discount rate, these sensitivity runs yield relatively similar 
SC-CO; distributions with mean values that differ by -18% to +11% (Extended Data Table 1) from 


our preferred SC-CO» estimate (Extended Data Fig. 1). 


The single largest contributor to the overall increase in the SC-CO; relative to the widely used 
DICE model is the use of a lower near-term discount rate, with updated damage functions being 
the second largest contributor. We disaggregate impacts of the changes to the near-term 
discount rate, the sectoral damage functions, and the remaining GIVE components (the RFF-SPs 
and FalR) in Table 1. We-start by running DICE-2016R, which uses none of our updated 
components and uses DICE's default discounting approach, yielding an SC-CO» estimate of $44/t- 
CO». Updating the climate modelling, the socioeconomic scenarios, and the discounting approach 
reflecting a 3% near-term discount rate but retaining the DICE-2016R damage function increases 
the mean SC-CO» by 34% to $59/t-CO;». Incorporating our sectoral damage functions in place of 
the DICE-2016R damage function further increases the estimate to 580/1-СО;, or a total increase 
of 8196. Finally, using a lower 296 near-term discount rate has the largest effect, increasing the 
mean SC-CO; estimate to this study's value of $185/t-CO2, a 32196 increase relative to $44/t-CO2, 


and a 3.6-fold increase relative to the widely cited US government value of $51/t-CO;. 


The four climate damage sectors represented in the model vary substantially in their respective 
contributions to the overall magnitude and uncertainty of the SC-CO» (Fig. 3). Temperature 
mortality impacts are the largest driver of the SC-CO», contributing a mean partial SC-CO2 
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(defined as the SC-CO2 estimated for an individual impact sector) of $90/t-CO2 ($39-165/t-CO;, 
5-95% range) to the $185/t-COz total using a near-term 2% discount rate. Agricultural impacts 
have a similar mean contribution of $84/t-CO», but greater uncertainty, with a 5-95% partial SC- 
СО; range spanning -523 to $263/t-CO2. This large range, which includes the potential for 
beneficial effects of higher temperatures and CO» concentrations on agriculture, arises due to 
compounding uncertainty in the relationship between CO», temperature, and crop yields and 
how these factors interact with the economic system to affect human welfare!*. We sample 
uncertain parameters for mortality and agriculture (see Methods), the damage sectors for which 


parameter uncertainty is quantified in the underlying studies. 


The relatively small contribution of sea-level rise, which includes both coastal damages and 
adaptation costs, to the total SC-CO2 (mean partial SC-CO» of $2/t-CO2, $0-4/t-CO», 5-95% range) 
is attributable in part to the inertia in the physical system connecting CO» emissions and sea-level 
rise and in part to the optimal regional adaptation response allowed by the Coastal Impact and 
Adaptation Model (CIAM) that we incorporate into GIVE?!. Such optimal, forward-looking 
adaptation responses can substantially reduce estimated coastal damages relative to a static 
scenario assuming no response to evolving coastal гіѕк52°°6. Future research could improve the 
characterization of plausible versus optimal coastal adaptation responses. The relatively slow 
pace of sea-level rise also causes the greatest damages to occur far in the future when 
discounting effects are strongest. Energy costs for residential and commercial buildings (based 
on Clarke et al. 2018“) also make a relatively small contribution to the overall SC-CO2 (mean 
partial SC-CO» of $9/t-CO», $4-15/t-CO», 5-95% range), due to increased energy demand from 
cooling being offset by decreased heating demand and future technological progress; these 


results аге broadly consistent with other recent empirical work”. 


We quantify impacts on four critical, globally significant damage sectors that are often considered 
to contribute the most to the SC-CO;!^? and for which studies exist that can be readily 
incorporated into SC-CO» estimation due to their global coverage, regional disaggregation, and 
monetization. A limitation of this study is that other categories of climate damages, including 
additional non-market damages other than human mortality, remain unaccounted for. The 


inclusion of additional damage sectors such as biodiversity??, labour productivity995!, conflict®, 
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and migration® in future work would further improve our estimates. Current evidence strongly 
suggests that including these sectors would raise the estimates of the SC-CO», although 
accounting for adaptation responses could potentially counteract some of that effect. Other 
costs of climate change including the loss of cultural heritage, particular ways of life, or valued 
ecosystems, may never be fully valued in economic terms but would also likely raise the SC-CO2 
beyond the estimates presented here. The addition of alternate studies covering the same 
sectors to incorporate additional independent lines of evidence is also a promising area for 
continued work to improve the SC-CO2. The modular structure of the Mimi.jl framework 
facilitates such addition of new damage sectors with ease, providing a flexible basis for future 


scientific improvement of the SC-CO;. 


While we approximate the effects of a rapid Antarctic ice sheet disintegration tipping point within 
the BRICK sea-level component, incorporating additional potential discontinuities in the climate 
system would further improve our SC-CO; estimates9^^. We expect that, in total, the future 
inclusion of additional damage sectors and tipping elements is likely to raise the estimates of the 
SC-CO2, and therefore that the estimates from the present study are likely best viewed as 
conservative. Similarly, accounting for different climate model structures, as the recent IPCC AR6 
report does in chapter 7^, would further strengthen the robustness of our SC-CO» estimates and 
their associated uncertainty levels. For example, that chapter (see Cross-Chapter Box 7.1, Table 
2 therein) shows the MAGICC climate model projects slightly higher temperature increases than 


the FalR model. 


The methods employed in this study reflect the culmination of several important advances: 
development of fully probabilistic very long-run socioeconomic inputs that natively incorporate 
uncertainty over future climate policy; incorporation of state-of-the-science representations of 
the climate system and sectoral damage functions; and an empirically calibrated discounting 
approach that accounts for uncertainty in future economic growth. These advances collectively 
allow for the full characterization of uncertainties, and their compounding interactions, 
throughout all steps of SC-CO; estimation, including sectoral market and nonmarket damages to 
human health. Their implementation on Mimi.jl5, an open-source, modular computational 


platform for assembling IAMs, improves the scientific basis and transparency of the resulting 
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estimates and is responsive to the NASEM near-term recommendations. The methodology also 


provides a straightforward means to calculate SC-CO2 results for other years and estimate the 


social cost of other greenhouse gases (e.g., CHa, N20, HFCs). Our higher SC-CO2 values, compared 


to estimates currently used in policy evaluation, substantially increase the estimated benefits of 


greenhouse gas mitigation, and thereby increase the expected net benefits of more stringent 


climate change policies. 
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404 Tables 


Mean SC- Incremental Share of 


CO; Change Total 
Row Scenario ($/t-CO;) ($/t-CO;) Change (%) 
a DICE-2016R $44 
b GIVE w/ DICE damage function, 396 near-term discount rate $59 $15 1196 
c GIVE, w/ sectoral damages, 3% near-term discount rate $80 $21 15% 
This study: 
d GIVE, w/ sectoral damages, 2% near-term discount rate $185 $105 74% 


405 Table 1 | Evolution of mean SC-CO; from DICE-2016R to this study. All SC-CO2 values are expressed in 
406 2020 US dollars per metric tonne of СО». (a) represents the SC-CO» using base DICE-2016R deterministic. 
407 The mean SC-CO; of $44/t-CO; is similar to the value previously estimated from IWG DICE-2010 of $46/t- 
408 СО» at a 396 discount rate, after converting to 2020$°°, (b) then retains the DICE-2016R damage function 
409 but otherwise deploys GIVE under discounting parameters of p = 0.8%, 7 = 1.57, which are consistent 
410 witha 3% near-term discount rate, (c) then replaces the DICE-2016R damage function with our sectoral 
411 damage functions, (d) then uses our preferred discounting parameters from this study of p = 0.2%, 7 = 
412 1.24, which are consistent with a 296 near-term discount rate. The final row represents the preferred 
413 mean value from this study. 
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Figure legends 

Fig. 1 | RFF-SP socioeconomic scenarios and the resulting climate system projections. a-c, Probabilistic 
socioeconomic projections for global population (a), per capita GDP growth rates (b), and carbon dioxide 
emission levels (c) from the RFF-SP scenarios. d-f, corresponding climate system projections that account 
for parametric uncertainty in FalR and BRICK for atmospheric carbon dioxide concentrations (d), global 
surface temperature changes relative to the 1850-1900 mean (e), and global mean sea-level changes 
relative to 1900 (f). In all panels, solid centre lines depict the median outcome, with darker shading 
spanning the 25-75% quantile range and lighter shading spanning the 5-95% quantile range. 


Fig. 2 | $С-СО» distributions vary with the choice of near-term discount rates. Distributions of the SC- 
СО» based on RFF-SP scenario samples, a stochastic, growth-linked discounting framework, uncertainty in 
the FalR climate and BRICK sea-level models, and uncertainty in climate damage parameters. Colours 
correspond to near-term average discount rates of 3.0% (blue), 2.5% (orange), 2.0% (red, our preferred 
specification), and 1.5% (teal). Dashed vertical lines highlight mean SC-CO» values. Box and whisker plots 
along the bottom of the figure depict each SC-CO2 distribution’s median (centre white line), 25-7596 
quantile range (box width), and 5-95% quantile range (coloured horizontal lines) values. All SC-CO2 values 
are expressed in 2020 US dollars per metric tonne of CO;. 


Fig. 3 | Partial SC-CO? estimates and uncertainty levels strongly differ across the four climate damage 
sectors. Box and whisker plots for the climate damage sectors included in the GIVE model, based on partial 
$С-СО» estimates for each sector. Figure depicts the median (centre white line), 25-75% quantile range 
(box width), and 5-95% quantile range (coloured horizontal lines) partial SC-CO? values. Black diamonds 
highlight each sector's mean partial SC-CO2, with the numeric value written directly above. АП SC-CO» 
values are expressed in 2020 US dollars per metric tonne of CO». 
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Methods 


Socioeconomic projections 

The RFF-SPs? used in this study were designed to address the requirements for socioeconomic 
projections posed by SC-CO2 estimation: (1) The roughly 300-year time-horizon required.to 
account for the vast majority of discounted future damages; (2) the need for geographically 
disaggregated estimates of GDP and population to support damages at a regional scale; (3) 
uncertainty accounting for expected future changes in both technology and policy (the SC-CO; is 
measured against the best estimate of future emissions, inclusive of future mitigation policies 
except the one under analysis); and (4) the interdependence of future population, GDP, and 


greenhouse gas emissions trajectories‘. 


The RFF-SPs address key shortcomings identified in the approach to. socioeconomic projections 
originally developed by the US IWG in 20106 and used consistently through the current US 
interim estimates?. Тһе IWG used five socioeconomic scenarios to 2100, drawn from the Energy 
Modeling Forum 22 modelling exercise?", опе of which represented future climate policy. The 
IWG scenarios were critiqued for not spanning the true uncertainty in GDP, population and 
emissions, nor reflecting the broader scenario literature overall?^95, The RFF-SPs used here 
improve on those scenarios by explicitly characterizing uncertainty in the demographic, economic 


and emissions projections. 


The multi-century time horizon required for the projections is long relative to the length of the 
historical record available to estimate country-level statistical models of population and 
economic growth. Accounting for uncertainty in future emissions over that time horizon requires 
assessing the potential for structural changes in technology and policies that are out of the range 
of historical experience. To address these challenges, the RFF-SPs were generated based upon a 


combination of statistical and expert-based approaches. 


We generated probabilistic, country-level population projections through 2300? by extending the 
fully probabilistic statistical approach used by the United Nations (UN) for its official population 
forecasts to 2100. We further incorporated feedback and improvements suggested by a panel of 


nine leading demographic experts convened to review preliminary results. 
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Our trajectories of country-level GDP per capita Кот 2018 to 2300 come from a multifactor 
Bayesian dynamic model, in which each country’s GDP per capita is based on a global frontier of 
developed economies and country-specific deviations from that frontier!?. We reweight the 
probabilities of the Bayesian model trajectories using results from the RFF Economic Growth 
Survey, a formal expert elicitation focused on quantifying uncertainty in long-run economic 


growth?. 


The resulting probabilistic socioeconomic trajectories represent an alternative to existing 
scenario-based approaches, such as those based on the Shared Socioeconomic Pathways 
narratives. Such scenarios do not typically come with associated probabilities, though there have 
been efforts to assign such probabilities to the SSPs a posteriori based on expert surveys?’. The 
use of non-probabilistic scenarios have been criticized in the literature for being overconfident 
and failing to reflect uncertainty®’. Indeed, multi-century socioeconomic projections are deeply 
uncertain, as illustrated by the wide 5-95% ranges that we consider (see Figure 1). The scenarios 
based on the SSP narratives and their commonly used extensions beyond 210087072 fail to span 


that uncertainty.? 


We also generate multi-century distributions of global CO2, CH4, and N20 emissions through RFF's 
Future Emissions Survey, which elicited experts in socioeconomic projections and climate policy?. 
Experts provided uncertainty ranges for future fossil fuel and process-related СО» emissions as 
well as changes in natural CO» stocks and negative-emissions technologies, incorporating their 
own uncertainty around future mitigation policy. They also quantified the sensitivity of emissions 
projections to future economic growth, thereby allowing for the development of a joint set of 
projections of emissions and economic growth. The experts additionally provided uncertainty 
ranges for trajectories of CH4 emissions, N20 emissions, and net СО; emissions from other 
sources of CO» emissions and sinks. 

Climate models 

(1) FAIR 

We represent the global climate system and carbon cycle dynamics using version 1.6.2 of the 
Finite Amplitude Impulse Response (FalR) тоде!.7377 FalR is an emissions-based simple climate 


model with a carbon cycle that depends on background warming levels and cumulative carbon 
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uptake by land and ocean sinks. This state-dependency enables FalR to replicate the equilibrium 
and impulse-response behaviours found in more sophisticated Earth system models, which is 
important for producing scientifically grounded SC-CO» estimates. These features are not found 
in the previous climate models used for SC-CO» calculations, which lack carbon cycle feedbacks 
and have been shown to respond too slowly to changes in radiative forcing". We run FalR with 
randomly sampled CO», CHa, and N20 emissions time-series from the RFF-SPs and represent other 
greenhouse gases and short-lived climate forcers using the SSP2-4.5 scenario", the scenario that 
most closely matches the median RFF-SP emissions trajectories. We account for climate model 
uncertainties by randomly sampling a calibrated 2,237-member ensemble of parameters that 
was produced using FalR as part of the IPCC AR6”4. See Supplemental Information section SI.2 for 


more detail on the FalR model. 


(2) BRICK 

We make probabilistic projections of regional changes in sea level using the Building blocks for 
Relevant Ice and Climate Knowledge (BRICK) model. BRICK represents individual contributions to 
sea level from the Greenland and Antarctic ice sheets, glaciers and small ice caps, thermal 
expansion, and land water storage and has been thoroughly described in prior studies’. BRICK 
downscales changes in global sea level to regional changes using maps of time-invariant scaling 
factors!?77, The Antarctic ice sheet model component also accounts for a potential tipping point 
where rapid ice sheet disintegration can occur when annual mean Antarctic surface temperatures 


cross an uncertain threshold!$. 


We closely follow past work and calibrate BRICK to the historic sea-level record over the period 
1850-2017 with a Bayesian framework!^!^/97», This calibration process uses observational 
constraints on global mean sea-level changes?? in addition to individual contributions from 
glaciers and small ice caps?! , the Greenland ice sheet®**?, the Antarctic ісе sheet?^, and trends in 
thermal ехрапѕіоп?°. It further statistically accounts for measurement error estimates provided 
with each observational time-series data set*ó. We select physically informed prior distributions 
for BRICK's uncertain parameters that are consistent with previous model calibration studies!?7, 
For the Antarctic ice sheet model component, we select prior distributions based on a 


paleoclimate calibration that uses independent sea-level data from 240,000 years before the 
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current ега to the present?®. We use our calibration framework to create a Markov chain of ten 
million representative samples from BRICK's joint posterior parameter distribution and assess 
convergence based on graphical diagnostics and Gelman-Rubin potential scale reduction factors 
that are less than 1.1878, We discard the first one million samples for the initial burn-in period 
and select a random subset of 10,000 samples from the remaining chain for our final sea-level 
parameter values. The distributions of the uncertain parameters in BRICK are shown in 
Supplemental Information Table 4. 

Damage functions 

(1) Sea-level rise 

The sea-level rise damage calculations are based on the work of Diaz?! which presents the Coastal 
Impacts and Adaptation Model (CIAM). CIAM is an optimization model that assesses the costs of 
various adaptation strategies against flooding damages and potential impacts from regional 
changes in sea level. It chooses the least-cost strategy for each of over 12,000 coastal segments 
across the globe in the Dynamic Interactive Vulnerability Assessment (DIVA) database?? after 
taking into account local physical and socioeconomic characteristics. CIAM's potential adaptation 
strategies are specified as a combination of (1) a choice on retreating inland from the coastline, 
protecting coastal communities and infrastructure, or remaining in place without taking any 
adaptive actions and (2) a choice on the degree of investment in coastal defence against several 
different storm surge return periods conditional on protection being decided on. The DIVA 
database provides generalized extreme value distributions that define these return periods for 


each individual segment. 


CIAM is a-deterministic model. All uncertainty in coastal damages is therefore the result of 
uncertain sea-level projections that arise due to GIVE's probabilistic emission scenarios and 


climate and sea-level model parametric uncertainties that we sample. 


(2) Building Energy Expenditures 

The energy demand damage function is based on the results of Clarke et al. (2018)”°, a study that 
used the Global Change Analysis Model (GCAM)???! to project how climate change affects 
regional building energy demand through 2100. GIVE's damage functions relate each degree of 


global temperature rise to a change in regional energy expenditures, expressed as a proportion 
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of that region’s GDP. We derive these damage functions using output data provided by the 
authors of Clarke et al?°. That output includes, for each of the 12 GCAM regions, the net change 
in regional energy expenditures as a proportion of regional GDP at various temperature levels 
(varying over both time and scenario). Clarke et al.2° note that this relationship is approximately 
linear in temperature. For each of the 12 GCAM regions, we fit a linear function to these 
datapoints by regressing the net change in energy expenditures as a proportion of GDP on global 
temperature rise relative to the preindustrial period. We assume the intercept is zero to ensure 
the resulting function yields no change in energy expenditures at zero temperature rise. This 
yields a coefficient for each region, denoted В? (see Supplemental Information Table 2 for these 
values). Energy damages for each country i located in region j are then calculated using the 


corresponding coefficient, as 


Change in Energy Expenditures as a Proportion of GDP;, = В? х (Temperature Rise),. 


We multiply this energy expenditure share by country-level GDP to generate damages in dollars. 


Clarke et al.?? did not feature any explicit consideration of uncertainty, so we do not include 
uncertainty in this damage function. Uncertainty in energy-related damages remain, however, 


due to GIVE's uncertain temperature projections and GDP trajectories. 


(3) Temperature-related mortality 

The mortality damage functions are based on the results of Cromar et al. (2022)??, who convened 
a panel of health experts to conduct a meta-analysis of peer-reviewed research studying the 
impacts of temperature on all-cause mortality risk, which includes human health risks related to 
a broad set of health outcomes including cardiovascular, respiratory, and infectious disease 
categories. The-meta-analysis combined studies to produce regionally disaggregated estimates 
of the effects on all-cause mortality of each degree of warming across a broad range of baseline 
temperatures, including both increased mortality risk at high temperatures and reduced risk at 
cooler temperatures. This produced, for each of 10 regions, a point estimate (and its standard 
error) representing the net change in all-cause mortality risk per degree Celsius of globally 


averaged surface temperatures (see Supplemental Information Table 1). 
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To reflect uncertainty in these estimates, we sample these parameters Bj for region j from a 
normal distribution centred on the point estimate and set the standard deviation equal to the 
reported standard error. We then compute temperature-induced excess deaths in country i in 
region j as 
(Temperature Induced Excess Deaths) ;, 
= Bj х (Temperature Rise), x (Baseline Mortality); +, (2) 
where we calculate baseline mortality as the regional population level times its baseline mortality 


rate from the RFF-SPs, 

(Baseline Mortality) ;, = Population;,, x (Baseline Mortality Rate); ç. (3) 
We monetize these excess deaths using the value of a statistical life (VSL) as follows: 
Monetized Excess Mortality; = VSL;, x (Temperature Induced Excess Deaths); +. (4) 


The baseline VSL value for 2020 for the United States (denoted У51285 0) is derived using EPA’s 
1990 Guidance value of $4.8 million and adjusted for income growth and inflation, resulting in a 
2020 U.S. VSL of $10.05 million in 20205 (U.S. EPA, 2010) (see data explainer notebook in 
replication code for the full derivation). We then base the VSL for country i in year t on the EPA’s 


baseline VSL for 2020, adjusted for country i’s GDP per capita in year t, as 


en | GDP per capita, ) (5) 


VSLi è= VSL x 
p U$2070 " V GDP per capitays 2020 
where = = 1 represents the income elasticity of the VSL. The primary function of = is to adjust 
the US VSL to other countries and at uncertain future income levels. We use a unit elasticity which 


is in line with the central tendency of values recommended in the literature for such cases?^?^, 


(4) Agriculture 

The agricultural damage function is based on Moore et al. (2017), which estimated damages in 
two steps using: (1) a meta-analysis of published studies of the effects of temperature, rainfall, 
and СО; on crop yields that builds on previous work by Challinor et al. (2014)? and Porter et al. 


(2014)37; and (2) a computable general equilibrium model to estimate the economic welfare 
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consequences of these yield shocks while accounting for trade patterns and supply and demand 


adjustments in agricultural markets across 16 regions. 


Moore et al. (2017) present their results in the form of damage functions that directly relate 
global mean surface temperature increase to welfare change in economic terms. Their study 
presents three different parameterizations of these damage functions to characterize 


uncertainty: a central, low, and high estimate. 


They estimated each of these three parameterizations for 1, 2, and 3 degrees Celsius of 
temperature increase, resulting in three piecewise linear damage functions for each region (see 
Supplemental Information Figure 1). To address uncertainty as part of our Monte Carlo sampling 
framework, we sampled a value from a triangular distribution with lower bound 0, mode 0.5, and 
upper bound 1 for each draw. Assigning the low, central, and high damage functions to each of 
these values respectively, the two nearest functions were linearly interpolated to produce the 
damage function for that draw, also interpolating linearly between the resultant 1-degree Celsius 
value and the origin since damages at zero temperature increase can be assumed to be zero. 
Importantly, this uncertainty sampling scheme preserves the covariance between regions arising 


through connections in the global trade network. 
Lastly, we incorporated their results into our model via the equation, 
YPCit 


YPCi1990 
ag share 


-є 
AgPctCost;, = а ( ) fic), 


where AgPctCost;, is the damage in the agricultural sector as a proportion of GDP in region i at 
time t; o; is the share of agriculture in GDP in 1990 in region i; є = 0.31 is the income elasticity 
of the agriculture share in СОР; and f; is the piecewise linear function for region i resulting from 


the steps described above. 


Discounting 
Our discounting approach directly follows from NASEM recommendations as developed by 
Newell, Pizer, and Prest??!, Given the long residence time of CO» in the atmosphere, the damages 


from СО; emitted today persist for centuries. These future damages must be converted to 
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present dollar equivalents using an appropriate discount rate. The climate economics literature 
typically uses Ramsey-style discounting that links the discount rate to future economic growth’. 
This linkage leads to the Ramsey-like equation for the discount rate over time, denoted r,: r, — 
p +n X gt, where p is the rate of pure time preference, 4; is the average rate of consumption 
growth from the year of the emissions pulse (described in the next section) to year t, and 7] X g, 
reflects the extent to which society discounts damages because future individuals are relatively 
wealthier. More specifically, у reflects how much the marginal value of consumption declines as 
consumption increases (a 196 increase in consumption corresponds with a n% decline in the 


marginal value of a dollar). 


We evaluate the stochastic discount rate for each realized path of uncertain consumption growth 
( = p + nge), explicitly and structurally modelling the uncertainty in discount rates that is often 
summarized by a declining term structure!??, This uncertainty in the discount rate leads to a 
stochastic discount factor (SD Е,) used to discount future marginal climate damages. The SDF, 
can also be written equivalently in terms of relative consumption levels?^10! as 


son = SNC)". © 
(1+ p)! \c2020 


We use this SDF, to discount marginal climate damages (M D+) to a present value. 


While the climate economics literature routinely uses а Ramsey-like approach to 


discounting3254.101-105 


, prior estimates by the US IWG disconnected discounting and future 
economic growth-by using a constant, deterministic discount rate. That approach implicitly 
assumes that.;] = 0, corresponding to no linkage between consumption growth and discounting 
as well as zero aversion to risk. Our approach re-establishes the Ramsey-like link between growth 
and discount rates. We use p and у values that were empirically calibrated? to be consistent with 
the RFF-SPs and evidence on the observed behaviour of interest rates*®. This procedure also 
produces near-term risk-free discount rates (defined as the average risk-free discount rate over 
the first decade of the time horizon) consistent with the desired values, such as those reported 


in Fig. 1. Our preferred SC-CO» estimate corresponds to a near-term 2% rate, which is consistent 


with real risk-free interest rates over the last 30 years, and uses p = 0.2% and у = 1.24331. The 
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(p,n) values corresponding to the alternative near-term rates of 1.5%, 2.5%, and 3% are 
PN 


(0.01%, 1.02), (0.5%, 1.42), and (0.8%, 1.57), respectively. 


The Ramsey-like form for the discount rate is a standard approach to value marginal impacts and 
account for their risk amid uncertainty in future payoffs and consumption levels in the discounted 
expected utility framework^??^, In that framework, the value of the л parameter reflects the 
degree of risk aversion as well as the inverse of the intertemporal elasticity of substitution. That 
framework is also used for benefit-cost analysis of policy and regulatory analysis under 
uncertainty, as it quantifies the risk premium associated with uncertainty and risk aversion in the 
valuation of a marginal emission of CO». While the Ramsey framework is widely used, other 
considerations for decision-making under uncertainty in context of climate change, such as the 
role of epistemic uncertainty and alternative preference structures including ambiguity aversion, 
have also been proposed?. We use the discounted expected utility framework because it is the 


most established and widely used framework for regulatory and policy analysis19^108, 


Estimating the SC-CO; 
We estimate the SC-CO2 in a three-step calculation process. In the first step, we run the GIVE 
model out to the year 2300 for two separate cases: a "baseline" case and a "perturbed" case that 
adds an extra 0.1 МЕС pulse of СО; emissions in the year 2020 and is otherwise identical. In the 
second step, we calculate marginal»climate damages in year t as the difference in modeled 
damages per tonne between the pulse and baseline runs as 

4 Ка 


MD, = > У (Damages with Pulse,g, — Baseline Damages), (7) 


d-1r-1 
where we aggregate over each of the four damage sectors d at their respective geographic 


resolutions (i.e., countries or regions) r. 


In the third and final step, we calculate the SC-CO: by discounting these marginal damages using 
the stochastic discount factors SDF, from equation (5) above and then aggregate them over time 


into a single present value 
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2300 
SC-CO2 = 2: SDF, х МО, . (8) 


t=2020 


For our preferred results, we calculate 10,000 unique SC-CO2 estimates. For each estimate, we 
sample the RFF-SP scenarios to account for uncertainties in global CO2, CHa, and N20 emission 
trajectories in addition to country-level population and GDP growth levels. We also sample 
parametric uncertainties in the FalR and BRICK models as well as the agricultural and 
temperature-related mortality damage functions (Extended Data Table 2). As described above, 
our preferred SC-CO2 estimate uses discounting parameters of р = 0.2% and у = 1.24 for a 


near-term rate of 2%. 


When we report partial SC-CO2 estimates for a given damage sector, we follow the estimation 
procedure outlined above, but only include the impacts from that individual sector when 
calculating marginal damages in equations (7) and (8). We normalize our estimates based on the 
emission pulse size and report all results throughout the paper in units of 2020 US dollars per 
metric tonne of CO». We use the implicit GDP price deflator from the U.S. Bureau of Economic 


Analysis to convert values to 2020 dollars. 


We typically summarize the distribution of our 10,000 SC-CO» estimates by its mean, i.e., 
E[SC-CO2], where the expectation operator is taken jointly over all uncertain parameters 
determining marginal damages (MD,) and the stochastic discount factor (SD F+). This calculation 
is consistent with economic theory for pricing investments and other actions with uncertain 
payoffs, and therefore properly accounts for the risk premium in the valuation of a marginal 


emission of CO» owing to the many compounding uncertainties we тоде! 9. 


Software 
All our results are computed using open-source software tools. We use the Julia programming 
language for the entire replication code of this paper!??. All models used in this study are 


implemented on the Mimi.jl computational platform for integrated assessment models?. 
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Data and Code Availability 
The replication code and data for this paper is available at 
https://doi.org/10.5281/zenodo.6932028, including instructions on how to rerun the entire 


analysis for this paper. 
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Extended data figure/table legends 


Extended Data Table 1 | Mean SC-CO> values (with 5*-95* quantile ranges), by damage 
function and discount rate ($/t-CO;). Our preferred estimates correspond to the GIVE sectoral 
damage functions at a 296 near-term discount rate, shown in bold. All results use the RFF-SP 
scenarios, a stochastic growth-linked discounting framework, and sample uncertain climate, sea 
level, and damage function parameters, including for DICE-2016R and Howard & Sterner 
damage functions. The DICE-2016R damage function is based on Nordhaus 2016 (see page 2 of 
Nordhaus 2016 Supplemental Information)??. The Howard & Sterner damage function is based 
on the base coefficient in their Table 2, specification (8). All SC-CO» values are expressed in 
2020 US dollars per metric tonne of СО». 


Extended Data Table 2 | Sources of SC-CO; uncertainty. The left column shows the inputs and 
components of the GIVE model that contribute to uncertainty in the SC-CO;. The right column 
briefly describes these uncertainties and their sources. 


Extended Data Fig. 1 | SC-CO; distributions are robust to different damage function 
specifications ($/t-CO2). Distributions of the SC-CO2 using the damage functions from GIVE 
(orange, our preferred specification), DICE-2016R?? (blue), and Howard & Sterner?? (red) for near- 
term discount rates of 1.596, 2.096, 2.596, and 3.096. All results use the RFF-SP scenarios, a 
stochastic growth-linked discounting framework, and sample uncertain climate, sea level, and 
damage function parameters, including for DICE-2016R and Howard & Sterner damage functions. 
The DICE-2016R damage function is based on Nordhaus 2016 (see page 2 of Nordhaus 2016 
Supplemental Information)?. The Howard-& Sterner damage function is based on the base 
coefficient in their Table 2, specification (8). АЛ SC-CO? values are expressed in 2020 US dollars 
per metric tonne of СО». 


Extended Data Fig. 2 | Discounted marginal damages by year, preferred 296 near-term discount 
rate case. Solid line represents mean discounted marginal damages for a one-tonne CO; 
emissions pulse in 2020, dotted line represents the median, with darker shading spanning the 
25-7596 quantile range and lighter shading spanning the 5-9596 quantile range. All SC-CO» values 
are expressed in 2020 US dollars per metric tonne of CO». 
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Extended Data Table 1 


Model Component Uncertainty Source 


Global CO», CHa, and N2O emission trajectories RFF-SPs? 


Country-level GDP growth rates RFF-SPs?.10 

Country-level population RFF-SPs? 

FalR climate-carbon cycle model 2,237-member constrained ensemble of the uncertain parameters (sampled with replacement) 
from IPCC AR6 report” 

BRICK sea-level model 10,000-member ensemble of the uncertain parameters derived from a Bayesian calibration 
framework!5:15 

Agriculture damage function Uncertain damage coefficient distributions based on Moore et al.18 

Temperature-related mortality damage function Uncertain damage coefficient distributions based on Cromar et al.'? 
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